Optimized production of large Bose Einstein Condensates 
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We suggest different simple schemes to efficiently load and evaporate a "dimple" crossed dipolar 
trap. The collisional processes between atoms which are trapped in a reservoir load in a non adiabatic 
way the dimple. The reservoir trap can be provided either by a dark SPOT Magneto Optical 
Trap, the (aberrated) laser beam itself or by a quadrupolar or quadratic magnetic trap. Optimal 
parameters for the dimple are derived from thermodynamical equations and from loading time, 
including possible inelastic and Majorana losses. We suggest to load at relatively high temperature 
a tight optical trap. Simple evaporative cooling equations, taking into account gravity, the possible 
occurrence of hydrodynamical regime, Feshbach resonance processes and three body recombination 
events are given. To have an efficient evaporation the elastic collisional rate (in s _1 ) is found to 
be on the order of the trapping frequency and lower than one hundred times the temperature in 
micro-Kelvin. Bose Einstein condensates with more than 10 7 atoms should be obtained in much 
less than one second starting from an usual MOT setup. 

PACS numbers: 03.75.Ss, 32.80.Pj, 34.10.+X, 42.50.Vk 
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The creation of degenerate atomic and molecular gases 
is one of the major achievements in the last decade in 
physics. Using a Magneto Optical Trap (MOT) laser 
cooled sample as the starting point, the largest Bose Ein- 
stein Condensate (BEC) ever created contains 5 x 10 7 
atoms Ql The hydrogen BEC contains 10 9 atoms 
but starts from a cryogenically cooled sample [4]. The 
techniques to form a degenerate BEC or Fermi gas have 
evolved but it is still challenging to reach quantum degen- 
eracy in less than one second with large atom numbers. 
Some cooling methods such as superposed laser cooling 
HH0, Doppler cooling on forbidden transitions Q or 
degenerate Raman sideband cooling using for instance 
a three-dimensional far-off-resonant lattice have reached 
high phas e-sp ace density than MOTs but with complex 
setups HEflEJIIl- 

Up to now, the only way to reach the quantum de- 
generacy regime is to load atoms into conservative traps 
and perform evaporative cooling. This requires com- 
plex setups, both in terms of vacuum and in terms of 
optical access. The most common traps are magnetic 
traps, in which BEC can be achieved after tens of sec- 
onds of forced evaporation. The use of atom chips pro- 
vides a much tighter trap confinement than usual mag- 
netic traps. Consequently evaporative cooling is much 
more efficient, leading to BEC in less than a second but 
with atom numbers on the order of 10 5 13]. Another 
common trap is the optical trap. The use of CO2 dipole 
traps (also called QUEST trap) has led to large samples 
of degenerate gases of 5 x 10 5 atoms in few seconds start- 
ing from a MOT HI [H E E3 • Finally combinations of 
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traps have been partially studied. For instance the pow- 
erful so called 'dimple' trick, has been demonstrated [lif : 
a small but tight trap is superimposed on a large atomic 
trap acting as a temperature reservoir. For a loading effi- 
ciency of nearly 20 % of the initial number of atoms, the 
temperature remains almost unchanged and a large gain 
in the atomic density and in phase space density can be 
obtained, allowing fast evaporative cooling. 

Several theoretical studies of evaporative cooling exist 
(see [H 

and references therein) but none of them 
simultaneously takes into account gravity, the hydrody- 
namical regime effects, the ability to tune the scattering 
length value a through a Feshbach resonance, the temper- 
ature dependence of the scattering cross section a, three 
body recombination processes and the trap shape evolu- 
tion. In this paper, we present a study of combined traps 
as well as a theoretical model which takes into account all 
these effects. This provides tools that help devising the 
best strategy to reach quantum degeneracy, with many 
atoms, in a short time, and with a simple experimental 
setup. We will see for example that a simple magnetic 
plus crossed dipole trap setup can lead in much less than 
one second to a large (> 10 7 atoms) BEC from a standard 
vapor cell MOT setup. 

In the first part of this paper, we study loading strate- 
gies for a dipole trap from different types of reservoir 
traps. We mainly study the case of a dimple dipole trap 
focussed on a large magnetic trap. Dark-SPontaneous- 
force Optical Trap (dark-SPOT) [12 or degraded laser 
beam could also be used as a reservoir to load atoms into 
a dimple optical trap. We also derive thermodynamical 
equations for the loading process. We discuss the dy- 
namics of the loading process of the dimple, as well as 
the effect of Majorana losses and of two and three body 
inelastic losses. We emphasize that loading diabatically 
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the dimple trap is (three times) less time consuming than 
trying to load it adiabatically, and that it leads more 
rapidly to high phase-space densities needed for efficient 
evaporative cooling. 

The last part of this paper is devoted to the derivation 
of simple differential equations for evaporative cooling. 
We include gravity, hydrodynamical regime, the possibil- 
ity of modifying a with a Feshbach resonance, the tem- 
perature dependence of the scattering cross section, three 
body recombination and possible trap shape evolution in 
our model. Strategies for reaching quantum degeneracy 
are discussed. 



EFFICIENT LOADING OF AN OPTICAL TRAP 

Our goal is to provide a setup as simple as possible to 
rapidly reach quantum degeneracy by evaporative cool- 
ing. The speed of evaporation is directly linked to the 
trap frequency. If we want to avoid to work with atom 
chips the solution to have a tight trap is to use an op- 
tical trap, which also allows a better optical access for 
further experiments on the BEC or on the degenerate 
Fermi gases. In this article, we mainly focus on the load- 
ing and evaporative processes in a (crossed) optical trap. 
As mentioned in the introduction one of the most power- 
ful techniques to load an optical trap is to superimpose 
it on a large reservoir of cold atoms and wait for the col- 
lisions to fill the "dimple" formed by the optical trap. 
Once the dimple is loaded the reservoir trap is removed 
and the evaporation process is achieved by lowering the 
laser intensity to reach the quantum degeneracy regime. 

Before we describe in detail the case where the atomic 
reservoir is a magnetic trap, we would like to briefly dis- 
cuss other possible reservoirs, such as a dark-SPOT, and 
an aberrated optical trap. To use a dark-SPOT as the 
reservoir of atoms is an attractive idea for two reasons: 
first the atomic density is higher than in magneto-optical 
traps 0,H2- Second, at least for alkaly atoms, atoms 
are in the lowest hyperfine state, in which inelastic colli- 
sions are less likely, and they barely see the MOT light. 
Pioneer works have been done in a non degenenerate case 
[Hill, but very recently the group of Paul Lett @] at 
NIST reached quantum degeneracy with sodium atoms 
after loading a crossed dipole trap from a dark-SPOT. 
We believe that these results are due to efficient dimple 
loading directly from the dark-SPOT. In the experiment 
no depumping light is used, because the dark SPOT alone 
is specially efficient for sodium j2|| |27} , but we believe 
this might be needed for other atoms. 

A crossed optical trap alone (i.e. after that the MOT 
has been turned off) can also play the role of reservoir. 
Indeed, the efficiency of the loading process in a CO2 
laser cross trap has been clarified by the work on Ytter- 
biums atoms [2^, |2^| indicating that the atoms, trapped 
by the strong restoring longitudinal force existing in a 



CO2 trap, become concentrated into the cross 'dimple' 
region by atom-atom collisions. The CO2 laser is very ef- 
ficient because it diverges more rapidly (for similar waist) 
than smaller wavelength laser beams (Nd:YAG for in- 
stance), forming a longitudinal trap additionally to the 
radial one. The crossed region acts as a dimple in the 
reservoir formed by the non crossed region. The main 
drawback of using a CO2 laser is that it requires spe- 
cial windows (such as ZnSe one) in the setup. From the 



waist propagation formula w(z) = 



-Mw) 2 ' 



where the minimum waist is u>o — (created by a 

/ focal length lens and a collimated laser beam of waist 
W) , we see that a Far Off Resonance Trap (FORT) laser 
(A w 1.06 fim) such as a Nd:YAG or an Yb fiber laser with 
M 2 = 10 create the same trapping potential shape as a 
TEM 00 C0 2 laser (A = 10.6 /zm) with M 2 = 1. The near 
and far field beam shape of a CO2 M 2 — 1 laser cannot 
be both perfectly matched by a near infrared M 2 = 10 
laser but this is probably a small effect and is beyond 
the scope of our article. We then suggest that a CO2 
laser could be replaced by a M 2 sa 10 infrared laser, 
such as a nearly gaussian shaped multimode (fiber) laser 
or a diode laser. The trapping potential depth can be 
matched by using appropriate laser power. With similar 
trapping potential, if keeping a negligible photon scatter- 
ing, the efficiency of the loading and of the evaporative 
cooling should then be the same in both laser setup but 
with the simplicity of using a laser light non absorbed by 
the glass cell. 

Another possible way is to use a time averaged optical 
trap in a similar way to the one described in reference 
[20j. For instance by rapidly modulating the frequency 
driving an Acousto-Optic Modulator and thus by sweep- 
ing the position of the first order diffracted beam in front 
of a focusing lens it should be possible to increase the 
trapping volume of the reservoir without changing sig- 
nificantly the waist size at the lens focus neither its lo- 
cation (which defines the dimple). This time-averaged 
trap can also be used to optimize the spatial intensity 
shape during the evaporative cooling process. Finally, it 
should be possible to use aberrated laser beams (e.g. by 
computer-addressed holograms) [U H3 • 

We now turn to the use of a magnetic trap as a reservoir 
to load an optical trap, which is the focus of our study. 
Magnetic traps are the largest available traps with more 
than 10 9 atoms at T w 100 /iK directly transferred from a 
MOT Superconducting magnets can even catch, 

but from a buffer gas cooled, 10 12 atoms at 600 mK |34| . 
Recently, in reference H3, a Cr BEC has been obtained 
by using such a transfer technique from a magnetic to a 
crossed optical trap. Unfortunatly, in this experiment the 
depth of the optical trap was very small and, even after 
a RF magnetic trap evaporative cooling step, it was (in 
temperature units) only twice the atomic temperature. 
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This led to very fast decay due to plain evaporation after 
removing the magnetic trap. 

The goal of this article is to provide an optical trap 
with good loading and good starting condition before 
starting evaporative cooling. For simplicity reason, we 
will treat the case where two orthogonal laser beams 
cross in a central region, which is assumed to form a 
radial isotropic gaussian trap potential U(r) = £/i asor = 
— C/oe _2r / w ° . The potential depth Uq is then propor- 
tional to the one beam laser intensity / = where 
P is the laser power. With these simplifications the 
trapping angular frequency u> verifies u> = . where 



m is the atom mass (uj — * " ■> in the real physi- 
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cal crossed dipole trap situation of identical but orthog- 
onal laser beams). The experimental parameters: P 
and wo are related by Uq — r/ksT oc P/w 2 , oc lu 2 Wq 
where r/ = Uo/ksT is a dimensionless parameter, ks 
is the Boltzmann constant and T the temperature of 
the trapped sample. We illustrate the full potential 
(quadrupolar trap plus optical trap) in figure^for an op- 
tical trap with Uq — fcs x 1720/jK and u — 2tt x 1050 Hz. 
We will illustrate most of our results using the cesium 
atom (/ = 3, rn/ = —3 which has a magnetic moment 
fi ps — 3/is/4 > 0, where [Ib < is the Bohr mag- 
neton) because it is known to accumulate several prob- 
lems (mainly due the high value of its scattering length 
a = — 3000 do) |2a| such as hydrodynamical regime [37j |. 
large two body relaxation rate (3^| and large three body 
relaxation rate |3^| . In the cesium case the optical poten- 
tial in figure ^ can be created by crossing two Nd:YAG 
laser beams of 100 Watt each focused on a wq — 100/xm 
waist. 



Thermodynamics in an optical trap 

In order to study the loading of the dimple trap we 
will first deduce thermodynamics quantities (such as the 
atom number N and the temperature T) from standard 
thermodynamics equations 40]. Here, we will safely 
neglect gravity (see figure . Knowing T, N and a 
given potential U, we deduce the phase space density 
D = -7^- = uqA 3 where uq is the peak atomic density 

and A = 



— i — = is the thermal de Broglie wavelength. 
The D(N, T, U) function is calculated from the one atom 

U(r)-U(0) 

partition function Z\ = A~ 3 J Q 47rr 2 e k B T dr (for 
an infinite potential depth). We then deduce the Hel- 
moltz free energy F — NksT(h\D — 1), the entropy 
S = -|f and the energy E(N, T,U) = F + TS of the 
sample. For a power law potential, U(r) = U'r 3 / S and 
Z x = A- 3 (!^) S ^™ (r is the Gamma function) lead- 
ing, in the harmonic case where U(r) — ^muj 2 r 2 , to D = 
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FIG. 1: Trapping potentials (in fiK) for the cesium / = 
3,m/ = —3 atom. Top: Magnetic (thick black solid line) 
and crossed optical (thin blue solid line) potentials and its 
quadratic approximation (thin magenta dashed line). The 
magnetic gradient is lOmT/cm. The one beam laser is a 
100 Watt Nd:YAG laser focus on a wo = 100 /im waist. Bot- 
tom: Full potential (thick black solid line), with gravity in- 
cluded (thin magenta dashed line). The full potential for the 
magnetically untrapped state / = 3, m/ = +3 without grav- 
ity is plotted in a thin blue solid line. 
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The atoms are assumed to be initially in a mag- 
netic trap with an isotropic trapping potential U = 

u 



magn j 



with initial number of atoms Ni, temperature 
Tj, phase space density D i: entropy Si and energy Ei — 
E(Ni,Ti, £7 ma gn)- The gaussian laser trap can be super- 
imposed, to form the final potential Uf = t/ magn + t7i ase r, 
either in an adiabatic way, i.e. with a slow change in 
laser power and the system evolves with constant entropy 
Sf — Si, or in a sudden way with constant energy. Up 
to now, most of the experimental and theoretical works 
were focused on the adiabatic process 0, 0> EH • How- 
ever, it is faster to suddenly apply a small dimple trap 
and the difference in the final temperature (see figure |2J) 
or atom number is marginal compared to the adiabatic 
loading process. In this paper, we choose to discuss the 
diabatic loading process. 

The initial atomic spatial density distribution in the 

t/ maEn (r) 

magnetic trap is rii(r) = (rio)ie 
energy Ei becomes E[ = Ei 



3 T , and the initial 
/ °° Anr 2 rii(r)U\ ascr {r)dr 
when the laser trap is suddenly added to the mag- 
netic trap. After thermalization the final energy Ef 
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of the atoms in the potential Ut should be the same 
as the initial one E[. The magnetic and laser po- 
tentials have not the same potential energy at the 
center (C/ magn (0) = t/i aser (oo) = 0) then we have 
E f = E(N, Tf, Uf) + NUf(O). This leads to equation 
E(N, T f , U f ) + NU f (0) = E' t (N, T h U f ) which is used to 
determine the final parameters of the sample after a sud- 
den transformation (the equation in the adiabatic case 
is S(N,T f ,Uf) = S(N,T h Ui)) such as the final tem- 
perature Tf from which we deduce the final peak den- 
sity (no)f- The number of transfered atoms Nf left 
in the dimple trap after removing the magnetic trap is 

U f (r)-U f (0) 

Nf = J^ d 47rr 2 (no)/e kBT f dr where the dimple ra- 
dius Td is defined by by Uf(rd) = 0. 

Our numerical results are given in figure [21 (see also 
the top part of figure EJ) where they are compared with 
the one dimensional quadratic + two dimensional box 
model developed in reference and that, after small 
calculations, lead in our case to the following formulas: 
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ridTi < and r\v — y where 
ZijA? is the effective volume of the 

1 /3 

magnetic trap (I = V e is the magnetic trap size) and 



W i th a = urn = 

V e = Ni/(n )i 



Vd = Nf/(no)f is the effective dimple volume. 
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found that Vd is typically two times ( 



We 
where 
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is the value calculated using the quadratic ap- 
proximation of the gaussian trap. The right side approx- 
imations in equations hold when the key parameter 
rjve~ A / T - f i=y rive nd is small compared to 1. As indicated 
by the figures [3 and E3 the analytical equations Q in- 
spired by the model described in reference are found 
to be qualitatively correct and might be used as a first 
guide toward the optimized strategy. 

Equations seem to indicate that a very small r\v ^5 
e~ Vd but deep (% ^> 1) dimple trap would lead to large 
atom number loaded in the dimple, and to a very high 
phase space density. On the contrary, we will show that it 
is almost impossible to reach degeneracy using only such 
a very deep and small dimple trap because of hydrody- 
namical regime, two or three body inelastic or Majorana 
losses. 



Dynamics of the optical trap loading 

In this section, we give a physical picture of the mecha- 
nisms for the diabatic loading of the dimple trap. The di- 
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FIG. 2: Final phase space density (top) and temperature (bot- 
tom) at thermal equilibrium after a diabatic loading of the 
optical crossed dimple gaussian trap loaded from a reservoir 
containing 10 9 atoms at temperature 150 fiK. We consider dif- 
ferent dimple waists wo but the power is adjusted to keep the 
dimple depth (77^ « 11.5) constant. The red thick dashed 
lines represent formula Q, the thin blue solid lines are their 
approximations for small rive Vd and the thick black solid lines 
result of our thermodynamical model, where losses are not in- 
cluded. The final temperature after an adiabatic loading of 
the dimple is also indicated with a thin magenta dashed line. 



abatic loading relics on thermalization, and our descrip- 
tion clarifies how collisions rapidly fill the trap. Our es- 
timate for the diabatic loading time £i oa d is on the order 
of a few times the collision time in the magnetic trap, 
comparable to the estimate one given in |42| . This time 
is smaller than the timescale for an adiabatic loading, for 
which the trap should be applied slowly, so that thermal 
equilibrium is maintained throughout the loading. In this 
paper we choose to study the diabatic loading which is 
faster than the adiabatic one, while allowing similar final 
temperatures (see figure |2J) and still large gains in phase 
space densities. 

The detailed physics of the loading process is complex 
and depend on three parameters that we will study one 
after the other: the probability to be transfered into the 
dimple (or to be ejected) during a collision, the atomic 
flux entering the dimple region (depending on the col- 
lisions and oscillations times in the reservoir trap) and 
the collisional rate in the dimple region. To simplify, we 
consider here that r\d > 4 which is favorable for the sub- 
sequent evaporative cooling step and we suppose in this 



FIG. 3: Three types of collision occurring between atoms to 
load a dimple from a reservoir. The vertical dashed lines 
limit the Nd atoms in the dimple region (optically trapped or 
not) from the others. First: collision between two magnet- 
ically but not optically trapped atoms (black hatched ones) 
with kinetic energy w (77^ + ^)kgTi near the center. There 
is a probability p t ~ 0.5 to transfer one atom in the optical 
trap (process 1). Second: between one un(-optically)trapped 
atom and one trapped atom (empty red ones). There is a 
probability p cv « 0.15 to eject (evaporate) the trapped atom 
(process 2'). Third: between two trapped atoms (filled blue 
ones) having kinetic energy on the order of E% = |fcsTi, 
this thermalisation process could rarely lead to evaporation 
(process 3). 



FIG. 4: Equipotential lines for the magnetic and optical po- 
tential described in figure Q The coordinates z and p — 
\J x 2 + y 2 are expressed in )im (a part with the unphysical 
p < region is also drawn by symmetry for clarity). The 
dimple radius and the magnetic trap size I are indicated. 
A typical atomic trajectory is shown with a collision occurring 
slightly before the oscillatory time t OS c. 



section that the dimple trap is small enough such that 
the reservoir is unaffected during the loading: Nf <C N 
and ksTf « ksTi for instance (see figure^- 

The first parameter (probability to be transfered or 
ejected) depends on the energy change occurring during 
a collision. For a Boltzmann distribution, the proba- 
bility that during a collision a given atom (with a typ- 
ical kinetic energy Ek ~ f^B? 1 ) acquires an energy 

greater than r/ksT can be estimated to be p ( J = 

f(v)/f(0) = 1 + 2e-"vW^ - Erf(^) where f( V ) = 

e' k 

frth^f V E' K e~ kBT dE' K and Erf is the error function 
\19L |43j . Following notations of figure |31 let us first 
study collisions between two magnetically but not op- 
tically trapped atoms: these atoms of the reservoir have 
a total energy larger than the depth of the dimple trap 
because they have a typical energy + |) ksTi while 
in the vicinity of r = 0. During one of these collisions 
one atom is transfered into the dimple (process 1 in fig- 
ure 0) if its energy becomes less than the trapping en- 
ergy rjdksTi, implying the other atom acquires an energy 
higher than (r]d + 3)ksTi which occurs with the probabil- 
ity ^ = p f fa±|l|g^ ) ~ p (1) ~ 0.5. If one atom from 

the reservoir collides with an optically trapped atom, the 
chance that one atom ends up trapped in the dimple is 
even higher (process 2 in figure |3J). In order to be con- 
servative we will assume pt ~ 0.5 during the full loading 



process. 

The reverse processes (2' and 3 in figure|3J), namely the 
ejection of an optically trapped atom after a collision, 
must also be taken into account. This probability p ev 
can be evaluated in the worst case (2') when the collision 
occurs with an atom from the reservoir, whose kinetic 
energy is ~ rjdkBTi. The dimple atom kinetic energy is 
~ ksTi, except during the negligible time of the very 
beginning of the loading when the dimple atoms are not 
thermalized. This collision is then similar (in the center 
of mass point of view) to a collision of two atoms with 
^dksTi/2 average kinetic energy inside a trap of rjdksTi 
depth. The collision has therefore a probability estimated 

t0 be P { vTksTih ) = p ^ ~ °' 3 t0 leave one atom ( thc 
originally trapped or not) out of the dimple. Thus, the 
ejection probability is p cv 0.15. As mentioned this 
value is evaluated in the worst case, we then assume in 
the following p cv <C 1 and neglect the marginal ejection 
process during the whole loading time. 

Once the collisional processes are known we should 
consider the atomic flux entering in the dimple region. 
We use a single notation Nd for all the atoms in the 
dimple region (optically trapped in the dimple or not, 
see figure |3J| . At the very beginning of the loading pro- 
cess the Nd ~ Nijt w riir\ value comes from atoms of 
the reservoir and, when the sample is thermalized the 

r 3 

Nd ~ Nf Ni-p value mainly comes from optically 
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trapped atoms. 

Two typical times have to be taken into account to un- 
derstand the dynamics of the loading. First, the "oscilla- 
tory" time tose, defined by mvf ~ ksTi where v% ~ l/t osc , 
after which an atom is almost back at the same position. 
This corresponds to the period in a harmonic trap. The 
second time is the average time interval between two col- 
lisions in the reservoir, t rn u rs — - — where a is the scat- 
tcring cross section and where the exact numerical factor 
depends on the potential trap geometry and will be dis- 
cussed in the section devoted to the evaporative cooling 
study. In absence of collisions, a fraction ~ r^/l 2 of the 
magnetically trapped atoms trajectories cross the small 
dimple radius as shown by the figure 0] (remember that 
p = \J x 1 + y 2 is in fact a two dimensional coordinate). 
To calculate the rate of atoms going through the dimple 
volume, one needs to take into account both this frac- 
tion, and the "ergodicity time" it is required to modify 
the trajectories. Indeed, once all the atoms whose tra- 
jectories go through the dimple are loaded into it, col- 
lisions in the reservoir are needed to modify the trajec- 
tories so that some more atoms have a trajectory going 
through the dimple. This time is of the order (exactly 
for a perfect isotropic harmonic trap) of the collisional 
time t co \\ in the magnetic trap |44| . For simplicity we as- 
sume that the magnetic trap is not in a hydrodynamical 
regime: t co u > t osc , which is slightly wrong in the po- 
tential we consider here (see figure ^1 where t osc = 9 ms 
and icon = 7 ms. With such assumption, the number of 
atoms passing through the dimple region during a time 
At > t coll is AiVP ass ~ N^^, 

Each atom passing in the dimple volume has a prob- 
ability ~ rd{ndcr) to collide during a single travel inside 
the dimple region and rid ~ is the density. Before col- 
liding with another atom from the reservoir, this atom 
passes through the dimple a maximum of j^a. times. 
Therefore, the probability pd to transfer such an atom 
is then pd ~ ptrd{ndO')^f aSL ■ Thus, the number of atoms 
transfered inside the dimple during At might be esti- 
mated to be ANd ~ AN^^pd leading to the differ- 
ential equation Nd ~ Nd wi? where the dot designs 
the time derivative. The loading process changes when 
the number of atoms inside the dimple region reaches 
K = n° d r d defined by ~ 1/ Pt = 2. Using 

t rn u w — = — w we found the relation nR = 2— n,-. 
When Nd > N9, the hydrodynamical regime occurs and 
any atom going through the dimple will collide with at 
least one atom of the dimple, the probability pd of trans- 
ferring an atom is then about one each time an atom 
enters the dimple region. The hydrodynamical regime is 
not favorable for evaporative cooling so we have advan- 
tage to choose a dimple such as Nf < N®. In this case 
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FIG. 5: Dimple loading time ti oa d (solid line), compared to the 
initial collisional time (dashed line) t co \\ (here 7 ms), (top) and 
number of collisions before the evaporation point p co u (bot- 
tom) which indicates when higher than unity the hydrody- 
namical regime (hatched area). The conditions in the loaded 
dimple trap are the same as in figure [5] Nd:YAG optical 
crossed gaussian trap (with constant depth) loaded from a 
lOmT/cm gradient magnetic quadrupole trap containing 10 9 
atoms at temperature 150 /iK. 

the loading time of the dimple is 

t load ^2W^ln(^)^2t coll ln(^) (2) 

This exponential time loading is eventually followed by 
a linear time loading (with pd = 1) if Nf > N®. The 
results for our dynamical loading model are illustrated 
in figure I5|78| . We could test the accuracy of our naive 
theory in two ways. First using the case of a dimple 
as large as the reservoir ~ I, which is not cover by 
the assumptions made to derive all the previous formula, 
but which is nevertheless useful to test if all the initial 
atoms are in fact found to be transfered into the dim- 
ple (remember that r\d > 4). In this case rif ~ ^T2 ni 
formula reads Nf ~ n f r d ~ 2nJ 3 w 2iVj which is only 
wrong by a factor 2. The second test is the fact that 
the hydrodynamical regime is found to be reached for 
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the same waist wq w 120 — 150 fim using our theory 
(see figure |B] where N d w TV/) or the exact parame- 
ter p co i] (see figure |SJ . It might also be useful to test 
here if the diabatic loading time ti oa d is smaller that 
the adiabatic one ii oa d,ad- If we suppose that the adi- 
abaticity criterion (for an harmonic trap) — <^ T e i, ^- 
is correct for a ratio value of 1/10 sa 1/12 (to simplify 
the final formula); this leads, in a non hydrodynamical 

Thus, 



regime where r o i = j^- < to, to Co 



latcoii ' 



Wad « 12*^1* (£). But B oc iV^/T^ leads to 



^-load.ad 



4i co ii In 



In 



2f 



load 



In m 



A 



reasonable estimated of the adiabatic loading time is then 

^load,ad ^ 3tl oa d- 

In conclusion, we found from figure [S] that a few reser- 
voir collisional times are enough to load the dimple in rea- 
sonable agreement with the estimate of |42j but we found 
a strong dependence on the dimple radius ra- We have 
advantage to choose a dimple radius rd such as the trans- 
fer atom number found by the thermodynamics equations 

Nf verifies Nf « N% ~ 2^ which is the maximum value 
(for an almost hydrodynamical regime in the reservoir) 
before reaching the hydrodynamical regime in the dim- 
ple. The final density in the dimple is then n/ rs ^72 ni 

reached after a loading time 2t co ii In f^^M . For an opti- 
cal dimple trap, the link between and the laser waist 
depends on rd and on the reservoir but r\ w 2iuq seems 

2 

to be a reasonable approximation leading to N d 4^ 
which is found (see figure [fJJ to be a very precise approx- 
imation. 



body collisions follows the equation N 2 ~ K 2 NdUd ~ 
-^2"d r d- Solving this equation leads to the final num- 
ber N 2 j of lost atoms due to two body collisions N 2 j ~ 



tco\\K 2 n 2 f r d . 



The lost fraction -J^ 

N f 



is then simply ~ t co \\K 2 rif , i-e. the final two body loss 
during one collisional time (which is the only relevant 
time during the dimple exponential loading). Similarly 
the three body atom losses (rate T 3 ~ L 3 n d ) lead to a 
lost fraction ~ |tcoll-^3^/ 5 i-e. roughly the final three 
body loss during one collisional time. 

One issue with a quadrupole trap is Majorana losses, 
due to spin flips occurring when an atom crosses the Ma- 
jorana sphere of radius rMaj ~ yj ^57 (typically on the 

order of one micron) where Vd is the velocity inside the 
dimple trap |49|. Since the spin flip occurs near a zero 
magnetic field the change in potential energy is marginal 
and can be neglected, as shown by comparing the thin 
solid blue line and the thick black solid line in the bottom 
inset of figure Q The density in the dimple trap is much 
higher than the one in the magnetic trap so spin flips 
will mainly occur from atoms already inside the dimple 
trap where Vd ~ -Jk^Tf /m. The spin flipped atoms are 
trapped in the optical trap. However, for these atoms, 
the combined magnetic plus dimple trap depth is finite 
(see figure ^) , which leads to evaporation at a rate r ev . 
This process is not necessarily a bad one because it might 
be seen as the beginning of an evaporative cooling inside 
the dimple. Collisions could also occur with non opti- 
cally trapped atoms with a maximum rate of p ev /t osc . 
Of course, another possible strategy, to combine small 
two body losses with negligible Majorana losses, is to fo- 
cus the crossed optical trap slightly off center in a region 
where the magnetic field is small but not zero. 



Losses during the optical trap loading 

We now describe how the previous results are modified 
in the presence of losses (inelastic losses and Majorana 
losses). 

We suggest to use a quadrupole magnetic trap as the 
reservoir because the trap volume is large and the os- 
cillation frequency is high. Then [/ magll ~ [iB'r where 
fi is the atomic magnetic moment and B' is the mag- 
netic field gradient. Furthermore, with near zero mag- 
netic field B at the center, where the optical trap takes 
place, a quadrupole trap is a good choice to limit the 
two body inelastic processes. Indeed, assuming that no 
resonance processes, such as shape or Feshbach resonance 
take place, the magnetic field dependence of the two body 
inelastic rate is, for the worst case, proportional to \f~B 
[H. HE \M, 0, IH For instance in the cesium case this 
rate is K 2 = 4n ~ 4 x 10" 11 B 2 (mT)r(^K)" a78 cm 3 s- 1 
where n is the average atomic density 38]. During 
the loading, the number N 2 of atoms lost due to two 



To give an estimate of the number of spin flipped 
atoms, we consider that no spin flip occurs through in- 
elastic collisions. This assumption holds here because 
for our quadrupolar trap case K 2 has a negligible ef- 
fect (see figure El . Following the previous reasoning con- 
cerning AiV^ ass we find that the spin flipped atom num- 



fMaj 



N 



r Maj At 



where 



ber during At is roughly N d ^ Wj 

tcoll.d ~ t co uni/nd is the collision time inside the dim- 
ple trap. This result modifies the Majorana rate cal- 
culation given in reference 49] by a factor As a 



consequence, l\ d 



N, 



' Maj 

-?n— 



■rid, which is exactly the 



n t rp coU 

equation corresponding to two-body losses, with K 2 re- 



placed by 



. Thus, the final spin flipped fraction is 



This strongly depends on the dim- 
ple radius, a large radius being better in term of losses. 

Taking into account two-body inelastic losses, as well 
as Majorana and three body atom losses we estimate the 



8 



number of lost atoms to be: 



Moss ~ Nf 



f 



Maj Tbf 



where we choose B « B'(a r )d when evaluating the K$ 
formula. Formula for the three body loss rate 1^3 , t co \\ 
or r ev are given in the next section. The number of 
loaded atoms Nf (calculated from the thermodynamical 
equation) might never been achieved if the number of 
atoms lost Moss is large, in this case we are not able to 
calculate the number of loaded atoms but we can estimate 
it to be Nf - AW 



Numerical results 

In this section, we discuss the strategy to efficiently 
load a dimple trap, and give numerical results in the 
case of cesium. For our numerical simulations, we as- 
sume as before a P — 100 W Nd:YAG laser focused on a 
wo = 100 jUm waist. We also assume a lOmT/cm gradi- 
ent magnetic quadrupole trap containing 10 9 atoms at 
a temperature of 150 fxK. This corresponds to a final 
rj = -TTT- ~ 8 (different from r ld = = 11.5) 



if _ k B Ti 

which is a convenient starting condition for evaporative 
cooling. 

Based on thermodynamical considerations, as illus- 
trated in figure |21 one may believe that a small volume 
dimple trap would insure a high initial phase space den- 
sity and a small heating and can even be compatible with 
large atom numbers. When three body losses, Majorana 
losses and the hydrodynamical regime are taken into ac- 
count the conclusion is quite different. This is especially 
true for cesium whose scattering length a rs — 3000 ao 
is large at zero magnetic field. It is essential to take 
into account Majorana losses, as well as two body and 
three body losses, which only depend on the scattering 
length a. For further evaporative cooling, one also needs 
to avoid the hydrodynamical regime, when the collisional 
probability p co \i ~ * coll,d is greater than one, i.e when en- 
ergetic atoms collide before being evaporated. The strat- 
egy is then to choose a dimple radius r d in order to load 
2 

Nf < 2N. t -f atoms in a dimple to avoid the hydrody- 
namical regime, while avoiding Majorana and inelastic 
losses, and then to start evaporative cooling. 

Numerical results for loading cesium atoms in a dimple 
are shown in figures and If from the pure thermo- 
dynamical point of view a small dimple trap looks to 
be the best choice, figure indicates that its loading is 
strongly affected by Majorana transitions and three body 
losses and figure [5] clearly indicates that the high density 
reached leads to the hydrodynamical regime and to such 
high three body recombination rate that it is problem- 
atic for subsequent evaporation. Finally the (diabatic) 
loading of the dimple is fast. We also conclude that in 
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FIG. 6: Number of transfered atoms and losses during the 
loading of the optical trap from the reservoir containing 
Ni — 10 9 atoms at temperature Ti = 150 fiK. Top: final num- 
ber of transfered atoms Nf resulting of our thermodynamical 
model where no losses are included. The red thick dashed 
lines represent formula fT). the thin blue solid lines are their 
approximations for small r\ve rid and the thick black solid lines. 
N%, which is a naive threshold atomic for the hydrodynamical 
regime (see text) , is also indicated with a thin magenta dashed 

w 2 

line as well as its approximation 4— a . Bottom: Lost fraction 
^° s ' of the loaded atoms due to three body relaxation (thin 
blue solid line), two-body inelastic collision (red dashed line) 
or Majorana transition (thick black solid line). To be visible 
we have multiplied the negligible two-body inelastic losses by 
a factor 100. 



order to avoid three body losses and the hydrodynamical 
regime, and to have large number of atoms Nf ~ 4i/7q/<7 
loaded, it is necessary to use a rather large waist. This 
has a cost in the final phase-space density, but, as we will 
see in the next section, evaporative cooling starting with 
such conditions is so fast that this remains the best strat- 
egy. In our case we conclude that a wo = 100 /im waist, 
2 x 100 W laser (crossed laser of 100 W each), is a good 
compromise. More than 10 8 atoms can be transfered in 
hundreds of millisecond, leading to a phase-space density 
of 1/700 before starting evaporation. 
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Our calculations show that the dimple loading due to 
collisions between atoms is a very efficient process which 
leads to very high phase-space densities compared to the 
ones in a MOT. In general, when % w 4 — 12, it is pos- 
sible to find a waist value for which 20 % of the atoms 
are loaded without significant heating, and still avoiding 
the hydrodynamical regime as well as the Majorana and 
the two or three body losses. Such large atom numbers 
loaded in the dimple are an ideal starting point for evap- 
orative cooling. After polarization of the atoms in their 
lowest state in energy, to avoid two-body inelastic losses, 
evaporative cooling is performed by lowering the opti- 
cal trap depth. This will be discussed in the following 
section. 



EVAPORATIVE COOLING 

In this section, we give a complete set of simple equa- 
tions for evaporative cooling, and we discuss the optimal 
evaporative cooling strategies. Our th eory ge neralizes 
the one developed in [^3, E3 (see also |UE3)j by in- 
cluding the effects of gravity, three-body recombination 
events, the possible dynamical modification of the scat- 
tering length by a Feshbach resonance, as well as a time- 
dependent shape of the dipole trap. We derive simple 
differential equations for the atom number N, the tem- 
perature T and the energy E of the sample trapped in a 
time dependent trap, that can be easily solved numeri- 
cally. 

We will assume an isotropic power law potential: 
U(r) = U'r 3 / S for r < r v (U(r v ) = r]k B T is the ini- 
tial trap depth). The optical trap case (rv — w/V%) 
is idealized by the quadratic potential (U(r) = |mw 2 r 2 
with S = 3/2). For a non isotropic potential, similar 
equations can be written, for instance in the harmonic 
case by taking the trapping angular frequency to as the 
geometric average of the three different axis angular fre- 
quencies [soj . 

We will focus on the high r\ regime (77 > 4) where evap- 
orative cooling is more efficient (for evaporative cooling 
at low trap depth see |43j). In this high n regime, we 
assume that all the corrective factors (see |5fj) such as 
Pi(rj) = 1 — T[l, rj\/T[l] are equal to one, where T is Euler's 
gamma function. This leads for instance to neglect the 
spilling terms in the evaporation equations. We assume a 
Boltzmann gas behavior which is a good approximation 
except very close to the quantum degeneracy regime. For 
Bos e ga ses, see [3 and for degenerate Fermi gases, 
see 0, ■ We will assume that quantum degener- 

acy is reached (for the real gas) when the phase space 
density D reaches 1 for the classical gas. 

We assume s-wave elastic scattering with an elastic col- 
lision rate r c i = riQ(jv r where v r = 4-i/^^ is the average 

Q 2 

atomic relative velocity, a — 1 °^ a ^ is the energy de- 




FIG. 7: Equipotential lines for an isotropic gaussian opti- 
cal trap including gravity (corresponding to two crossed 5 W 
Nd:YAG lasers focused on w ~ 100 /im). Azimuthal z and 
polar coordinates p = \J x 2 + y 2 in fim have been used (a 
part with the unphysical p < region is also drawn by 
symmetry for clarity). A "trajectory" with azimuthal angle 
9 is indicated crossing the evaporation (dashed) circle with 
rj'(8)kBT minimal evaporation energy. The evaporation ener- 
gies rjkBT,r)' (6)kBT, (77 + a g )kBT are indicated respectively 
for 0, 8, it/2 angles. 

pendent scattering cross section with k = To take 

into account an average of the energy dependent scatter- 
ing in the potential, the authors of reference [2(| modify 
the cross section formula. The equations (53)-(54) of 
reference [2 (it are complex. Here, we modify them by re- 
placing k 2 by k 2 v — k 2 1 j in the scattering cross section 
formula. We found that such an approximation is valid 
to an accuracy of 20 percent, and is more physically in- 
tuitive, as it explicitly introduces an average momentum 
k cv in the new a = a yi formula[79j For a three di- 

mensional evaporation, the evaporation rate is given by 
rj? v = r c ie _? '(r/ — (| + <5))/V2 and the average energy 
taken by the evaporated atoms is E£ v ~ (?/ + l)A^fcsT. 
These formula are 15% accurate compared to the ones 
containing the Pi (7/) corrective factors • 

Gravity effect 

Taking into account gravity, the full potential is 
U g (r, z) = U (r) — mgz. The evaporation efficiency is sub- 
stantially reduced by the gravitational field when a g = 
> 1 [HElia. For 5 < 3 the minimum of the po- 
tential U~ is no more in r = and we have to redefine the 
trapping depth by r^ksT = U g (r — rjj , z — ru) — U~ . To 
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be used in the evaporative cooling equation a g has to be 
written as an 77 function. For instance, in the harmonic 
trap case: 



mg I g 



mgru 



k B T k B T 



2-qk B T 



(3) 



l4r]k B T / 1 



1 / V2mw g 

2 V vk B T 



4U^ 

rnwx 



For the crossed dipole trap case, wq = y/2ru, oj 

and Uq (x P/wq oc oj 2 Wq expressions are still valid (but 
Uq = U(rjj) is no more the trap depth rjk B T except in 
the zero gravity case) . 

During a typical optical evaporation process the laser 
power is reduced and the effect of gravity increases. 
Gravity tilts the potential, and modifies the evaporation 
rate, as is illustrated in figure [7| For the sake of simplic- 
ity, here, we only consider the case of isotropic traps, and 
we only describe "central" (r = 0) collisions. After such 
a collision, a given atomic trajectory oscillates between 
z < and z > so atoms will escape toward the z > 
region where the potential is the lowest. We label the tra- 
jectory with the azimuthal angle 9 (such as z — rcosQ), 
between and 7f/2. Due to isotropy the density probabil- 
ity of the azimuthal angle between 9 and 9 + d9 is sin 9d9. 
The evaporate atom energy is mainly between the thresh- 
old energy ri'(9)k B T = U g (r = r u ,z — ruCOs9) — U~ 
and (rj'(6) + l)k B T [l^. Therefore, and for simplicity 
reason, we assume that the 9 angle is determined exactly 
by the evaporation energy. Then we choose to average 
the escape rate to r ov = JJ 

age energy to E cv = J Q -ft—E^ 
r]'{9) = i] + a g {\ — cos 9) are: 



t/2 tVW smQdO and the aver- 
, which with 



sin 



mean free path (na) 1 where n = n$2 5 is the average 
atomic density. In the harmonic trap case £ ps \pKo r 
mil, leading to the approximate formula p co u ps -rfj . 
This is similar to the crossover criterium between the 
hydrodynamical and the Knudsen regime given by |64| . 
We then approximately take into account the collisional 
chain by multiplying the atom and the energy losses by 

an empirical smoothing function f(p co i\) = (l +Pcoll) 
which smoothly tends to zero in the high hydrodynamical 
regime. We believe that our results are quite insensitive 
to the exact function which is chosen to describe hydro- 
dynamical effects because we tried to remain in situation 
where these effects are negligible. 



Three body recombination 



Three body recombination (TBR) leads to losses at 
a rate r 3 = L^n 2 — L^n 2 ^^ 8 . The detailed expres- 
sion is complex and depends on the molecular poten- 
tial. Here, we will use the upper bound expression: 



L 3 w 225^ 1+Q a 1(fca)4 which is an approximate formula 
based on the experimental results of reference |39l| and 
on theoretical works j6f| |(5(| |U E3 ■ This expression is 
valid if no resonance processes are present (such as Efi- 
mov states occurring for a < at very low temperature 
It leads (we use k ev instead of k in L3 formula) to 
a simple result: L3 < 3.9— a 2 and therefore 



T 3 < 0.15 x 3 3/2 



-5 



fiT c i 
k B T 



el- 



(4) 



Nk B T 



1 



- 1 



- 1 



We could easily check that the limit of no gravity a„ 







is correct. These formula are probably not accurate for 
strong gravity effects, but they will be used to test if 
the gravity is negligible in a given evaporative cooling 
scheme. 



Hydrodynamical effects 

Evaporation is efficient only in a non hydrodynami- 
cal regime where the atoms do not collide before reach- 
ing the evaporation point 61] . The probability p co \\ 
for colliding before the evaporation point can be esti- 
mated by dividing the typical size of the cloud £ by the 



Moreover, TBR leads to "anti-evaporation" (the loss of 
the coldest atoms). A spatial averaging of the poten- 
tial energy U(r) with the TBR rate T 3 leads to an en- 
ergy heating of ^5k B T per event. TBR formes molecules 
in their highest vibrational state, with binding energy 
krtTh — 2h „-> (formula valid for large scattering length 



3ma- 



value). We will not take into account complex possible 
ro-vibrational relaxation. But, we have to take into ac- 
count the fact that the atomic (or molecular) products 
of TBR are trapped in the sample, which is therefore 
heated by T h per event, if k B T h < r}k B T . 39]. This 
heating does not occur if k B Th > r\k B T because these 
particles are ejected from the trap, except if the hy- 
drodynamical regime is reached. To take into account 
this possibility, we consider that TBR leads to a heat- 
ing term of k B T^ multiplied by a smoothing function 

/tbr = l-/(Pcoli) (l - / (ff|#)) which is 1 in all cases 
except when k B Th > r]k B T and in a non hydrodynamical 
regime. 
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Evaporative cooling equations 

We consider a generic loss term T\ oss = Tin + Tbg where 
Tb g is the background collisional rate and Ti n = K^n is 
the inelastic loss rate. Ti n can be zero in an optical trap 
where the atoms can be trapped in their true ground 
state. The energy E of the atoms is modified by the 
change of the trap shape through the change of poten- 
tial energy E pot — SNksT = E 3 ^+s occurring at a rate 

Tpot = jj7 0, 113, . In an optical trap E can change 
by Nks times the recoil temperature T reco ii due to pho- 
ton absorption at a rate ri ascr oc r/ksT. Finally, the 
evaporative cooling equations to be solved are 

n = - [r ev /(p coll ) + r loss + r 3 ] n (5) 
e = — pioss + r 3 ] e - E cv T cv f(p co n) + i? P otr po t + 

,2 
3' 



N 



Tlascr -^recoil 



We stress that these equations take into account all 
relevant physical phenomena for evaporative cooling: 
temperature dependent cross-section, one-body (such as 
background collisions) and two-body losses, TBR, hydro- 
dynamical regime, gravity, arbitrary modification of the 
trap depth and shape, time-dependent scattering length. 
All terms in this set of equation can be written as a func- 
tion of N, T, a and U (i.e. r\ and u in the harmonic 
case) . Considering that — = ^ + ^ these equations lead 
to differential equations for N and T which only depend 
on a and U . Simple analytical solutions exist in simple 
cases, when there is no three body process, no gravity, a 
is not temperature dependent... In our case, these equa- 
tions need to be solved numerically, and we use a simple 
Mathcmatica program. These programs, as the one de- 
scribing the dimple loading, are available upon request 
to the corresponding author. 



Evaporative cooling strategy 



-p -^recoil -p 
1 laser ^ J- 3 



\ST + T h f TBR 
T 



(6) 



Tpot does not appear in this equation because changes 
in the potential shape do not affect the phase space den- 
sity. This formula is valid (as well as equations (J5J) when 

the adiabatic conditions <C T i and jp <C 55- (in the 
harmonic case) |52l |72( are verified. We have numeri- 
cally verified that all the results presented in this article 
follow these adiabatic conditions within 90%. Equation 
(JHJ) indicates that D/D is maximized for 77 w 4.1 + S, 
high r e i just before the hydrodynamical regime (optimal 
value r e i ~ 3w in the harmonic case) and negligible loss 

rate. We have as typical value ^ r+ ^ h/TBR » § + r/ « 10. 
Therefore, the TBR heating rate is ~ 10r 3 . We then 
see, from equation ©, that in order to have an efficient 
evaporation the TBR heating rate ~ 10r 3 , the loss rate 
Tioss, and the laser absorption rate ri a serIr eco ii/T must 
not be higher than r ov . Using equation (J3J and (in the 
harmonic case) T^ 6 = 0.0035r o i we found that the TBR 
term is negligible if Tc^s" 1 ) < 300T(/xK). An efficient 
evaporation should then verify for 77 « 6 



T^s- 1 ) w oj < 300T(/_tK) 



(7) 



For rj = 10 this becomes r c i(s _1 ) <C 15T(/iK) which is 
much harder to achieve. This equation already indicates 
that, in contrary to the common intuition, the choice of 
a high initial temperature is not a bad choice because 
it allows to have a tight trap with fast evaporation due 
to high collisional rate without reaching the three-body 
collisional regime. 

A very crude estimate of the minimum time iBEC 
needed to reach the degeneracy starting from an initial 
phase space density Di can be estimated from equation 
© with no losses and r c i w lu kept constant during the 
evaporation process leading to £bec ~ ~ 'nn?'" 1 f° r ?7 ~ 6 



or tr 



10 can be used to 



1 bec » -oMti for 1 ~ 10 - V 
improve the final number and rj = 6 is a faster strategy 

to reach BEC. 



The optimization of evaporative cooling is usually 
achieved by maximizing at each time 7 = (7 > 3 



for good experimental conditions |35|,|41j). But, if there 
is no inelastic, TBR or background gas processes, 7 is 
maximized for an infinite 77 value which leads to in- 
finite evaporation time [toL Flf . A better parameter, 
at least for the speed of the evaporation process, is 

D 

D N V2 

(D/D 



(| + <5)y + <5r pot . More fancy parameters 



such as 



N/N 



can also be used to optimize at each time 



the evaporative cooling process. ^ becomes, in the zero 
gravity case, 



Numerical results 

Before we detail the conclusions of our study, in this 
paragraph, we compare our numerical simulation to the 
experiment described in reference |l4| and modeled in 
reference We choose the parameters of the table 

1 of reference |5lj : rn is here the rubidium mass, a — 
100 d , A^o = 6.7 x 10 5 initial atom number, T = 38/j,K 
initial temperature, ujq — 2tt x 1500 Hz initial angular 
frequency, rr = 6 s and rj = 10. Formula J3Jl lead to an 
effective waist of w$ = 40 /im slightly different than the 
experimental waist of < 50 /urn. 

Our results are shown in figure [S] The hydrodynam- 
ical regime is never reached. When both gravity and 
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FIG. 8: Evolution of the phase space density during the evap- 
oration at 77 = 10 and constant laser waist (see text). The 
four curves represents respectively, from the bottom one to 
the top one, full numerical resolution of the equations 10, 
resolution with no gravity (g = 0.01 m 2 /s), without taking 
into account the F3 term, and with no gravity neither three 
body loss terms. 

TBR are neglected, we recover the theoretical results of 
[51j . However, when we take those effects into account, 
we predict that the choice of parameters of reference [Hl| 
would not lead to BEC. We found that it is neverthe- 
less possible to reach BEC, in agreement with the exper- 
imental observation [l^j . starting the evaporation with 
the experimentally measured parameters No = 2 x 10 6 
and To = 75 \xK. Our numerical treatment is therefore 
able to reproduce the experimental observations, which 
theories without gravity or TBR cannot do as accurately. 

Evaporation of the dimple optical trap 

We now turn to our numerical results and discuss the 
best strategy to reach BEC in the case of cesium. We 
start evaporation in conditions close to the one deduced 
from our dimple loading theory with a 100 W Nd:YAG 
laser focused on 100 ^m: N = 10 8 , T = 200 ^K, rj = 9. 
Such a dimple trap has a heating rate of ri asC r x T rcco ii w 
11 s _1 x 0.2 [iK at full power. This heating rate will be 
found to be negligible. All parameters T],(J,a,g can be 
experimentally adjusted, and our theory takes this into 
account. In particular, gravity g can be effectively sup- 
pressed by adding a vertical levitatin g m agnetic field gra- 
dient on the order of tens of mT/cm |7j|. This levitating 
magnetic field gradient also produces a parabolic anti- 
trapping potential, with frequencies on the order of a few 
Hz, which is usually a small effect that we will neglect 

As illustrated in figure 03 a strategy with constant 
waist ends up with gravity issues, due to quasi one dimen- 
sional evaporation. We illustrate different evaporative 



cooling strategies in figure using a levitating magnetic 
field but a constant waist to avoid the gravity problem 
and in figureEDusing a time changing waist. We assumed 
losses due to collisions with hot atoms from the back- 
ground gas, leading to a maximum lifetime of 3s = jy— • 
We choose a scattering length value of 100 ao because it 
is found to be close to the optimal value and, despite the 
mass difference, it can mimic the rubidium case. Even in 
the levitating magnetic field case, BEC is not reached if 
77 = 9 is kept constant. By rapidly reducing r\ to a value 
r\ = 6, three body recombination and hydrodynamical 
regime are avoided, and BEC is reached in one second. 
This picture points out the effects of TBR and hydrody- 
namical regime, which our theory can predict and help 
avoiding. 

In the left part of figure we show how to reach de- 
generacy in less than one second by dynamically changing 
the waist of the dimple. The strategy chosen here is to 
keep the trapping frequency u> constant (see equation [SJ, 
In this case, the condensated atom number is very small, 
because of hydrodynamical regime or three body losses. 
However, if a is dynamically modified, as illustrated by 
the right curves in figure^J| BEC can be reached in much 
less than one second, with large atom numbers. A Fesh- 
bach resonance is used to have a high collisional rate r c i. 
However, this rate is kept lower than 3a; to avoid the 
hydrodynamical regime and lower than 40 x T(/zK) to 
avoid TBR (the factor 40 depends on the rj value chosen, 
here rj « 9). We reach degeneracy in a 0.35 s evaporation 
ramp with N — 2 x 10 final atom number at a temper- 
ature of 14 fiK. During this time the waist is divided by 
a factor 3. 

To find 'the' optimized strategy to reach BEC by dy- 
namically changing all P, w, a or gravity parameters is 
beyond the scope of this paper. However, our studies 
show that keeping r e ] w u> w lOOT(^tK) for r\ w 8 pro- 
duces even bigger degenerate samples in less than 100 ms, 
by maximizing the collision rate, while avoiding the hy- 
drodynamical regime. 

CONCLUSION 

In this paper, we have studied the sudden superimposi- 
tion of a dimple trap on a large volume trap. We propose 
three schemes to experimentally realize this goal: a dark- 
SPOT, a (aberrated) dipole trap, or a magnetic trap as 
the atom reservoir. The dynamics of the loading leads 
to an exponential grow of the atoms number in time and 
the diabatic loading is three times faster than the adia- 
batic one. The optimal dimple trap is the tightest possi- 
ble trap before reaching the hydrodynamical or the three 
body recombination regime allowing an almost tempera- 
ture invariant diabatic loading, i.e. with a maximum of 
20% of the atoms loaded. We have shown that the ana- 
lytical equations can be used to find the right range 
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FIG. 9: Efficiency of the evaporative cooling from numeri- 
cal resolution of the equations |0 for constant waist, without 
gravity [g = 0.01 m 2 /s), a = 100 a , No = 10 s , T = 200 /jK. 
Left: 77 = 9, Right: r\ decreases in a 0.2 s exponential decay 
from 77 = 9 towards an asymptotic value of 77 = 6. Top: evo- 
lution of the atom number N (solid line) and the temperature 
T (dashed line). Middle: evolution of the phase space density 
D. Bottom: hydrodynamical p co ii (thin blue solid line), three 
body r3/r cv (thick black solid line) and gravity a g (dashed 
red thick line) dimensionless parameters. 



FIG. 10: Efficiency of the evaporative cooling from numeri- 
cal resolution of the equations for a waist zoom, chosen 
here to keep the trapping frequency lo constant, No = 10 8 , 
T = 200 fj,K, r) = 9. Left: a = 100 a and right: a = 30 a 
decreases in a 0.2 s exponential decay towards an asymptotic 
value of 10 ao- Top: evolution of the atom number N (solid 
line) and the temperature T (dashed line) . Middle: evolution 
of the phase space density D. Bottom: hydrodynamical p co n 
(thin blue solid line), three body r3/r cv (thick black solid 
line) and gravity a g (dashed red thick line) dimensionless pa- 
rameters. 



for the dimple volume and radius. The focus of our study 
is a magnetic reservoir and a laser crossed dipole trap for 
the dimple. The magnetic setup is very simple because 
only two independent coils can provide magnetic field for 
the MOT, the quadrupole reservoir, the levitation and 
even provide the field for a Feshbach resonance. A linear 
quadrupole trap can be chosen to avoid two body decay 
but the size of the dimple trap has to be large enough to 
avoid the effect of Majorana losses which have been stud- 
ied. Spin flip process can also be avoided if the crossed 
optical trap is slightly off axes, where the magnetic field 
is not zero, or by using a quadratic magnetic trap. In this 

case of a dimple created by crossed lasers the laser power 

2 

and waist Wq have to be chosen to load ~ 4— 11 (where 
a is the scattering cross-section) atoms from a magnetic 
trap in a quasi hydrodynamical regime. 

We then studied in detail teh evaporation inside the 
dimple trap. For the first time, to our knowledge, 
we have derived simple evaporation equations including 
two-body, three-body inelastic collisions, hydrodynami- 
cal regime, effect of gravity, possible dynamical modifica- 
tion of both the temperature-dependent scattering cross- 
section and of the trap parameters. These equations (J5J 



can be extended to anisotropic potentials, low dimen- 
sional gases or to gases with anisotropic electrostatic in- 
teraction Q • We discussed the best evaporative cooling 
strategy for rapidly reaching quantum degeneracy with 
large atom numbers. We suggest to use a tight trap, and 
an elastic collisional rate r e i as high as possible but not 
higher than 300T(/jK) (here 77 sa 6, but the factor 300 
depends on the 77 value chosen) to avoid the three body 
collisional losses and not higher than u) (in the harmonic 
trap case) to avoid the hydrodynamical regime. The 77 
parameter could be chosen between 6, for fast evapora- 
tion, and more than 10 for large final number of atoms. 
We described how a levitating magnetic field gradient 
could be used as well as dynamical modifications of the 
scattering length and the trap shape (laser waist for in- 
stance) to lead to a very fast evaporation and large final 
atom numbers. Our main conclusion, which also holds 
for Fermi or Bose mixtures [T^, is that, using 

a tight trap with high frequency to and high collisional 
rate r o i but low scattering length and high temperature, 
large (> 10 7 atoms) degenerate samples could then be 
achieved in much less than one second starting from a 
standard MOT. In our group an experimental activity is 
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under way to test these schemes. 
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